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Abstract. The large time behavior of nonnegative solutions to the 
Q , reaction-diffusion equation dfU = — (— A)"/^u — u^, {a G (0,2], p > 1) 

' posed on and supplemented with an integrable initial condition is 

studied. We show that the anomalous diffusion term determines the 
large time asymptotics for p > 1 + a/N, while nonlinear effects win if 
p < 1 + a/N. 



0^ 



< 

^ • 1. Introduction 



> 



We study the behavior, as t ^ oo, of solutions to the following ini- 
tial value problem for the reaction-diffusion equation with the anomalous 
diffusion 



: (1) dtU = -A"u + Xu"", xeR^,t>0, 

OS; (2) uix,0) = uoix), 

CN ■ where the pseudo-differential operator A" = (— A)"/^ with < a < 2 is 

defined by the Fourier transformation: A"m(^) = Moreover, we 

O ■ assume that A G {—1, 1} and p > 1. 

Nonlinear evolution problems involving fractional Laplacian describing 
^ ■ the anomalous diffusion (or a-stable Levy diffusion) have been extensively 

c3 ■ studied in the mathematical and physical literature (see [21 [HI |5] for ref- 

erences). One of possible ways to understand the interaction between the 
anomalous diffusion operator (given by A"" or, more generally, by the Levy 
diffusion operator) and the nonlinearity in the equation under consideration 
is the study of the large time asymptotics of solutions to such equations. 
Our goal is to contribute to this theory and our results can be summarized 
as follows. For A = — 1 in equation ([1]), nonnegative solutions to the Cauchy 
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problem exist globally in time. Hence, we study the decay properties of the 
mass M(t) = f^^f u{x,t) dx of the solutions u = u{x,t) to problem ([I])-®. 
We prove that hm^^oo = Moo > for p > 1 + a/N (cf. Theorem 1, 
below), while M{t) tends to zero ast— >ooifl<p< 1 + a/N (cf. Theorem 
2). As a by-product of our analysis, we show the blow-up of all nonnega- 
tive solutions to ([I])-(l2l) with A = 1 in the case of the critical nonlinearity 
exponent p = 1 + a/N (see Theorem 3, below). 

The idea which allows to express the competition between diffusive and 
nonlinear terms in an evolution equation by studying the large time be- 
havior of the space integral of a solution was already introduced by Ben- 
Artzi & Koch [Ij who considered the viscous Hamilton- Jacobi equation 
Ut = Au— \ (see also Pinsky [16]). An analogous result for the equation 
Ut = Au + |Vm|^ (with the growing- in-time mass of solutions) was proved 
by Laurengot & Souplet [13] . Such questions concerning the asymptotic be- 
havior of solutions to the Hamilton- Jacobi equation with the Levy diffusion 
operator were answered in [11]. 

In the case of the classical reaction-diffusion equation {i.e. equation ([Tj) 
with a = 2), for p < 1 + 2/N, Fujita [6] proved the nonexistence of nonneg- 
ative global-in-time solution for any nontrivial initial condition. On other 
hand, if p > 1 + 2/N, global solutions do exist for any sufficiently small non- 
negative initial data. The proof of a blow-up of all nonnegative solutions in 
the critical case p = 1 + 2/N was completed in [HI [171 [12]. Analogous blow- 
up results for problem ([I])-([2l) with the fractional Laplacian (and with the 
critical exponent p = 1 + a/N for the existence/nonexistence of solutions) 
are contained e.g. in fTfl ITj 18113] . 



2. Statement of results 

In all theorems below, we always assume that u = u{x,t) is the non- 
negative (possibly weak) solution of problem ([I])- ([2]) corresponding to the 
nonnegative initial datum uq G L^(R^). Let uq ^ 0, for simplicity of the 
exposition. We refer the reader to [5] for several results on the existence, the 
uniqueness and the regularity of solutions to ([I])- ([2]) as well as for the proof 
of the maximum principle (which assures that the solution is nonnegative 
if the corresponding initial datum is so). 

First, we deal with the equation ([T]) containing the absorbing nonlinearity 
(A = — 1) and we study the decay of the "mass" 

(3) M{t) = / u{x,t)dx= / uo{x) dx — / / vF{x,s)dxds. 

Remark. In order to obtain equality ([3]), it suffices to integrate equation 
([T]) with respect to x and t. Another method which leads to ([3]) and which 



NONLINEAR EQUATION WITH FRACTIONAL LAPLACIAN 



3 



requires weaker regularity assumptions on a solution consists in integrating 
with respect to x the integral formulation of problem ([I])-® (see (fT2|) . 
below) and using the Fubini theorem. □ 

Since we limit ourselves to nonnegative solutions, the function M{t) 
defined in ([3]) is nonnegative and non-increasing. Hence, the limit Moo = 
limt^oo M{t) exists and we answer the question whether it is equal to zero 
or not. 

In our first theorem, the diffusion phenomena determine the large time 
asymptotics of solutions to ([I])-©. 

Theorem 1. Assume that u = u{x,t) is a nonnegative nontrivial solution 
of (HD-dl) with A = -1 andp> 1 + a/N. Then limt^oo M{t) = > 0. 
Moreover, for all q & [I, oo) 

(4) t^(^"^)||M(t) - MooPa(t)||g ^ ast ^oo, 

where the function Pa{x,t) denotes the fundamental solution of the linear 
equation Ut + A°m = (cf. equation Q below). 

In the remaining range of p, the mass M{t) converges to zero and this 
phenomena can be interpreted as the domination of nonlinear effects in 
the large time asymptotic of solutions to ([Il)-(l2]). Note here that the mass 
M{t) = /jgjv u{x,t) dx of every solution to linear equation ut + A'^m = is 
constant in time. 

Theorem 2. Assume that u = u{x,t) is a nonnegative solution of problem 
(QD-dD with A = -1 andl <p<l + a/N. Then limf^oo M{t) = 0. 

Let us emphasize that the proof of Theorem 2 is based on the so- 
called the rescaled test function method which was used by Mitidieri & 
Pokhozhaev {of. e.g. [HI [15] and the references therein) to prove the nonex- 
istence of solutions to nonlinear elliptic and parabolic equations. 

As the by-product of our analysis, we can also contribute to the theory 
on the blow-up of solutions to ([I])- ([2]) with A = +1. Recall that the method 
of the rescaled test function (which we also apply here) was use in [71 [8] 
to show the blow-up of all positive solutions to ([l])-([2]) with A = 1 and 
p < 1 + a/N. Here, we complete that result by the simple proof of the 
blow-up in the critical case p = 1 + a/N. 

Theorem 3. // A = 1, a G (0,2] and p = 1 + a/N, then any nonnegative 
nonzero solution ci/([2])-(0) blows up in a finite time. 
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3. Proofs of Theorems 1, 2, and 3 

Note first that any (sufficiently regular) nonnegative solution to ([I])- ([2]) 
satisfies 

(5) ^— u{x,t)dx= / Uo{x) dx + X / / vF{x,s)dxds. 

Hence, for A = — 1 and Uq G L^(M^), we immediately obtain 

(6) u E L°°([0, oo), L\R^)) n LP{R^ x (0, oo)). 

Proof of Theorem 1. First, we recall that the fundamental solution Pq, = 
Pq,(x, t) of the linear equation dtu + A^m = can be written via the Fourier 
transform as follows 

(7) P„(x, t) = r^/-P„(xt-V«, 1) = ^_L_ e"-«-*l«l" dC 
It is well-known that for each a G (0, 2], this function satisfies 

(8) P„(l) G L°°(M^) nLi(M^), Paix,t)>0, I P^{x,t)dx = l, 



for all X G and t > 0. Hence, using the Young inequality for the convo- 
lution and the self-similar form of Pa, we have 

(9) \\Pa{t)*uo\\p < Cr^^'-'M/-llu,\\,, 

(10) l|VP,(t)||p = c'r^^^-i/^')/"-^/", 

(11) \\Pa(t)*Uo\\p < \\Uo\\p, 

for all p G [1, oo] and t > 0. 

In the next step, using the following well-known integral representation 
of solutions to (^-(^ 

(12) u{t) = Pa{t) * uo - [ Pa{t- s)*uP{s)ds, 

Jo 

we immediately obtain the estimate < u{x,t) < Pa{x,t) * Uo{x). Hence, 
by and we get 



(13) < min{Ct-^(^'-i)/"||^io||?;lko|0 =//(t,P,«,Mo). 

Now, for fixed e G (0, 1], we consider the solution = u'^{x,t) of ([1])- 
([2]) with the initial condition euo{x). The comparison principle implies that 
< M^(x, t) < u{x, t) for every x G and t > 0. Hence, it suffices to show 
that for small e > 0, which will be determined later, we have 



= lim / M=(x,t) dx > 0. 
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Note first tlie using equality ([5]) in the case of the solution , we obtain 

(14) M^ = e\ [ Uo{x)dx-- [ [ {u'{x,t)y dx dt\ . 

Urn ^ Jo Jk^ J 

Now, we apply (fT3l) with u replaced by u'^. Observe that the function H 
defined in (IT^ satisfies H{t,p,a,euo) = ePH(t,p,a,Uo). Hence 

- / / {u^{x,t)Y dxdt < - H{t,p,a,euo) dt 
£ Jo Jm £ Jo 



/■oo 

g-i-p / H{t,p,a,uo) dt. 
Jo 



It is follows immediately from the definition of the function H that the 
integral on the right-hand side is convergent for p > 1 + a/N. Consequently, 

- / {u^{x,t)Y dx dt ^ as e \ 0, 
^ Jo Jrn 

and the constant given by f|T^ is positive for sufficiently small e > 0. 

From now on, the proof of the asymptotic relation (jlj) is standard, hence, 
we shall be brief in details. First we recall that for every uo G L^(M^) we 
have 

(15) lim ||P„(t)* Mo- MP„(t) 111 = 0, 

t— >oo 

where M = f^j^ uo{x) dx. This is the immediate consequence of the Tay- 
lor argument combined with an approximation argument. Details of this 
reasoning can be found in ^ Lemma 3.3]. 

Now, to complete the proof of Theorem 1, we adopt the reasoning from 
|13j . It follows from the integral equation ( |T2l) and inequality (3.7) with 
p = 1 that 

\\u{t) - P^{t - to) * u{to)\\i < [ \\u{s)\\Pds for all t > to > 0. 

J to 

Hence, using the triangle inequality we infer 

\\u{t) - M^P^{t)\\i < [ IHsWpds 

J to 

(16) + WPcit - to) * u{to) - M(to)Pa(t - to) 111 

+ ||M(to)(Pa(t-to)-Pa(t))||l 

+ ||P„(t)||i|M(to)-Moo|. 

Applying first ([8]) and (fT5|) with uo = u(to), and next passing to the limit 
as t — s> oo on the right-hand side of fll6p . we obtain 

/>oo 

limsup||u(t) -Mo,P,(t)||i < / \\uis)\\lds+\M{to)-M^\. 

t— »oo Jto 
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By letting to go to +00 and using ([6]) we conclude that 

(17) \\u{t) - MooPait)\\i ^ as t^oo. 

In order to obtain the asymptotic term for p > 1, observe that by the 
integral equation (fT2l) and estimate ([9]), for each m G [1, 00], we have 

(18) ||n(t)||„ < \\Pait)*uo\\m < Ct-^(^-i/™)/"||^io||i. 
Hence, for every q G [1,"^), using the Holder inequality, we obtain 

\\u{t) - M^P^{t)\\g < \\u{t) - M^P^{t)\\\-' {\Ht)\t + WM^PaimL) 

< ct-^(^-i/^)/"||iz(t) - MooP«(^)||}-^ 

with S = (1 — l/g)/(l — 1/m). Finally, applying f|T7j) we complete the proof 
of Theorem 1. □ 

Proof of Theorem 2. Let us define the function ip{x, t) = {(pi{x)Y' {(p2{t)Y 
where 

" ^^"^ Vii^) = i^(^], V2{t)=ij(^], R>0, 



p-1 ^ ' ' " \BRJ ' ' " \R» 

and ip is a smooth non-increasing function on [0, 00) such that 



1 if < r < 1, 

if r > 2. 



The constant > in the definition of (pi is fixed and will be chosen later. 
In fact, it plays some role in the critical case p = 1 + a/N only while in the 
subcritical case p < 1 + a/N we simply put B = 1. In the following, we 
denote by fli and ^2 the supports of ipi and (p2, respectively: 

^]l = {x G : \x\ < 2BR] , ^2 = {t G [0, 00) : t < 2/?°} . 

Now, we multiply equation ([T]) by ^{x^t) and integrate with respect to 
X and t to obtain 

UQ{x)(p{x,{i) dx — I I u^{x,t)(p{x,t) dxdt 

u{x, t)ip2itY^''iMx)Y dxdt 

(19) - / / u{x,t)<pi{xYdt<p2{tY dxdt 

Jn2 Jqi 

< i / u{x,t)<p2{tYvi{xY'^^"Vi{x) dxdt 
u{x, t)ipi{xY^2{tY'^'^t^2{t) dxdt. 

^2 -Til 
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In ( |T9l) . we have used the inequahty A"y)f < ^ip^]~^A°'ipl, (see [H Prop. 2.3] 
and [lOj, Prop. 3.3] for its proof) which is vahd for all a G (0, 2], £ > 1, and 
any sufficiently regular, nonnegative, decaying at infinity function yji. 

Hence, by the e- Young inequality ah < ea^ + C{e)h^~^ (note that 1/p + 
l/(£ - 1) = 1) with e > 0, we deduce from ^ 

(20) 

Mo(x)(y9(x, 0) — (1 + 2^5) / / u^{x,t)(p{x,t) dxdt 



< c{e)e 



[ [ ipiifiilA'^ipif-^ dxdt+ [ [ ^{^2\dtV2t^ dxdt] . 



Recall now that the functions (pi and v?2 depend on i? > 0. Hence chang- 
ing the variables ^ = R~^x and r = R'^'t, we easily obtain from ( |20il the 
following estimate 
(21) 

uo{x)ip{x,0)dx-{l+2£e) [ [ uP{x,t)ip{x,t) dxdt < CR^+''-''^^-^\ 

where the constant C on the right hand side of ( l2TI) is independent of R. 
Note that N + a — a{£ — l) < if and only if p < 1 + a/N. Now, we consider 
two cases. 

For p < 1 + a/N, we have N + a — a{i — 1) < 0. Hence, computing 
the limit oo in fl2T]) and using the Lebesgue dominated convergence 

theorem, we obtain 



Moo = / uoix)dx- / / uP{x,t)dxdt <2ie / / dxdt. 
Jri^ Jo Jri^ Jo Jri^ 

Since u G Lp{R^ x (0, oo)) (cf. (ED) and since e>0 can be chosen arbitrary 
small, we immediately obtain that M^o = 0. 

In the critical case p = 1 + a/A^, we estimate first term on the right hand 
side of inequality (fT9|) using again by the e- Young inequality and the second 
term by the Holder inequality (with p = p/{p — 1) = i — I) as follows 

uq{x)(p{x,0) dx — / / uP(p{x,t) dxdt 
Ui Jn2 Jqi 



< ie / uP{x,t)dxdt 

Jq.2 J^l 

(22) +C{e) [ I ^f{t)^f-^'^\x)\h.^^i{x)fdxdt 

Jq.2 J^l 

+i{ I [ uP{x,t)dxdt) 

[ /fix)^t'^'it) \dtMt)f dxdt) 
n2 Jui / 
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Here, = {t e [0, oo) : R"" <t< 2R"} is the support of dt^2- Note that 

{x , t) dx dt ^ as R—>-oo, 



because u G Lp{R^ x [0, oo)) (cf. 

Now, introducing the new variables ^ = {BR)~^x, t = R'^^t and recall- 
ing that p = 1 + a/iV, we rewrite (j22l) as follows 

(23) 

uq{x)lp{x,Q) dx — I I uP{x,t)(p{x,t) dxdt — ei / / vF{x,t)dxdt 

i/p 



<CiB^^p([ [ uP{x,t)dxdt] +C2C{£)B 



where the constants Ci,C2 are independent of R, B, and of e. Passing in 
(123|) to the limit as i? — ^ +oo and using the Lebesgue dominated convergence 
theorem we get 



oo p roo 



Uo{x) dx — I I vP{x,t) dxdt — el I I vF{x,t)dxdt 
(24) jRiv Jo Jk^ Jo Jm^ 

< C2C{e)B-'^. 

Finally, computing the limit i? ^ oo in ^IM we infer that M^o = beacuse 
e > can be arbitrarily small. This complete the proof of Theorem 2. □ 



Proof of Theorem 3. The proof proceeds by contradiction. Let u be 
a non-negative non-trivial solution of ([Il)-(l2]) with A = 1. Take the test 
function the same as in the proof of Theorem 2. Repeating the estimations 
which lead to fl23ll. we obtain 



Mo(x)(/9(x, 0) (ix + / / vF{x,t)ip{x,t) dxdt 
^1 Jni 

— e£ / uP{x,t)dxdt 
Jn2 Jni 

<CiB^/p([ [ uP{x,t)dxdtY +C2C{e)B-''. 

Now, we chose e = l/(2£) in fl2Sl) and we pass to the following limits: 
first -R — s> oo, next B ^ oo. Using the Lebesgue dominated convergence 
theorem, we obtain 

1 f°° f 

uo{x) dx + - / uP{x, t) dxdt < 0. 
^ Jo Jr'^ 

Hence, u{x,t) = which contradicts our assumption imposed on u. □ 
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